Fine tuning of quantum operations performed via Raman transitions 
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A scheme for fine tuning of quantum operations to improve their performance is proposed. A 
quantum system in A configuration with two-photon Raman transitions is considered without adia- 
batic elimination of the excited (intermediate) state. Conditional dynamics of the system is studied 
with focus on improving fidelity of quantum operations. In particular, the it pulse and -n/2 pulse 
quantum operations are considered. The dressed states for the atom-field system, with an atom 
driven on one transition by a classical field and on the other by a quantum cavity field, are found. 
A discrete set of detunings is given for which high fidelity of desired states is achieved. Analytical 
solutions for the quantum state amplitudes are found in the first order perturbation theory with 
respect to the cavity damping rate k and the spontaneous emission rate 7. Numerical solutions for 
higher values of k and 7 indicate a stabilizing role of spontaneous emission in the n and n/2 pulse 
quantum operations. The idea can also be applied for excitation pulses of different shapes. 

PACS numbers: 03.67.Lx, 42.50.Ct, 42.50.Dv 
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I. INTRODUCTION 



An atomic A system with Zeeman sublevels of the 
ground state, which are sufficiently long lived to store 
qubits, plays a very important role in quantum compu- 
tations. Superpositions of such states can store quan- 
tum information for sufficiently long time to make var- 
ious quantum algorithms feasible Therefore, re- 
searchers frequentlyconsider such states in their propos- 
al, |a la , [7 1 |8| and use them in their experi- 
ToL I 111. H2, HH, HH . Access to a good quan- 
tum memory is necessary, but not sufficient, to realize 
quantum computation. It is also important to be able to 
perform quantum operations on such coded qubits with 
high fidelity differing from unity by 10~ 5 or less (TBI. Il6|. 

Although there is no direct, strong coupling between 
the Zeeman sublevels, one can efficiently manipulate pop- 
ulations of the states using the two-photon Raman tran- 
sition involving an auxiliary level, for example, state |2) 
in the A configuration, as shown in Fig. 1. This method 
of qubit encoding and manipulation has many different 
implementations. There are single atoms or ions mod- 
eled by different level configuration systems: a three-level 
A system [l2l . Il3j . a four-level system with three levels 
in A configuration, and one additional long-lived level, 
which does not participate in the transition (TtI [Lsj . a 
six-level double A system, which consists of two A sys- 
tems behaving exactly in parallel 2, 6, 19]. There are also 
solid-state implementations: quantum dots modeled by 
the A-type [20, [21( three-level system and superconduct- 
ing quantum interference devices modeled by the A-type 
three-level system [H H, H H HI . In all these sys- 
tems it is easy to store qubits and it is easy to perform 
single qubit gates just by turning the lasers on and off, to 
drive transitions for a proper period of time. Moreover, 



A-type systems are also perfect to realize the multiple 
qubit gates, as for example the crucial for quantum com- 
putation controlled NOT gate [H [13, [28[ , because one can 
place two or more such systems in a microcavity and cou- 
ple one transition in each A-type structure to the cavity 
field mode. Then, qubit interactions are mediated by the 
cavity field mode. Another way to realize two qubit gates 
is to perform the joint detection of photons leaking out 
of two separate cavities with trapped A-type systems [f| . 

Since there are many advantages of devices composed 
of A-type systems and cavities, such systems are very 
popular elements of various quantum information pro- 
cessing schemes. However, there is one important draw- 
back of such schemes — the population transferred to the 
intermediate level diminishes the fidelity of quantum op- 
erations performed with them. To avoid the destructive 
effect of populating the intermediate level, one can decide 
to work in the decoherence-free space @, 0, H3, HU, US H3] 
or, alternatively, one can assume that the cavity mode 
and the laser fields are far detuned from their respective 
transitions, as to make the population of the intermedi- 
ate state negligible. The latter approach is referred to as 
adiabatic elimination. In realistic situations with not ex- 
tremely large values of detunings, the population of the 
intermediate level is small but noticeable enough to sig- 
nificantly deteriorate the quality of quantum operations. 




FIG. 1: A-type Raman transition: a key ingredient of different 
quantum devices. 
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In this paper, we propose a scheme for fine tuning 
of the quantum operations performed via Raman tran- 
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sitions in the atomic A system driven on one transition 
by a classical field, and on the other by a quantum cavity 
field. We discuss the conditional evolution of the A sys- 
tem in a situation when the adiabatic elimination is not 
made, and we show that it is possible to take advantage of 
the fact that the intermediate state population oscillates 
rapidly, periodically approaching zero, which can be used 
to improve the quality of quantum gates based on the Ra- 
man transitions in the A system. It turns out that it is 
possible, by setting appropriate values for the detuning, 
to make the population of the intermediate state negli- 
gible at the time the quantum operation is completed, 
and in this way we are able to increase significantly the 
fidelity of the operation. Such fine tuning of the quan- 
tum operations works well even for small detunings and, 
therefore, it can prove useful when one wants to perform 
quantum operations with high fidelity. We have found 
a discrete set of detunings for which perfect operations 
are possible if there is no damping. We have also found 
approximate analytical solutions describing the evolution 
of the A system including both the cavity damping and 
spontaneous emission. Numerical results for higher val- 
ues of the cavity decay rate and spontaneous emission 
rate show, somewhat unexpectedly, that the spontaneous 
emission can play a stabilizing role improving the result 
of quantum operation when the cavity decay rate be- 
comes sufficiently large. We shortly address the issue of 
fine tuning for the nonrectangular pulses. 



II. MODEL 

Let us consider two basic operations, which can be 
performed on qubit-encoding states of the A-type level 
structures: the n pulse operation and the tt/2 pulse oper- 
ation. These two basic operations can be used to realize a 
number of important quantum information tasks as, for 
example, generation of maximally entangled states 
quantum information transfer [H, l3l| , or the performing 
of controlled two qubit gates @ . We consider an atom in 
A configuration trapped in a cavity. The long-lived states 
|0) and |1) of the atom are coupled via the intermediate 
level 1 2) (see Fig.[T]). The |0) <-> |2) transition is coupled 
to the cavity mode with a frequency w cav and coupling 
strength g. The second transition is driven by a classical 
laser field with the coupling strength Q. The frequency 
of the laser field is wl- Both the classical laser field and 
the quantized cavity mode are detuned from the corre- 
sponding transition frequencies by A = (E^— E±)/h— ojl- 
The population transfer between long-lived states |0) and 
|1) take place only when both transitions |0) «-> |2) and 
|1) <-> 1 2) are driven. Therefore, the qubit is safe if the 
laser is turned off. If we want to perform the tt pulse 
operation or the tt/2 pulse operation then all we need to 
do is turning the laser on for a proper period of time. 

The evolution of the A-type quantum systems is deter- 
mined by the effective non-Hermitian Hamiltonian (we 



set ft, = 1 here and in the following): 
H = (A — 17)022 + (^cr 2 i + gaa 2 o + H.c.) — ina)a , (1) 

where 7 is the spontaneous emission rate from the atomic 
state 1 2), and k is the cavity decay rate. In expression ((T|) 
we also introduce the flip operators <xy = where 
i, j — 0, 1, 2. The effective Hamiltonian ([I} describes con- 
ditional evolution of the atom-field system and will be 
used in our further calculations. 

III. EXACT SOLUTIONS 

As a first step, we find the solutions for the Schrodinger 
equation when both 7 and k are zeros. This assumption 
allows us to obtain exact solutions. The Hamiltonian ((T|) 
takes in this case the form 



H = Aa-22 + (fioai + 9 a(J 20 + H.c.) . 



(2) 



and it can be easily diagonalized in the basis being the 
product states of the atomic states and the cavity field 
photon states. To give the expressions a more compact 
form, we denote by \jn) — \j) ® \n) a state of the system 
consisting of the atomic state \j) and the cavity field 
with n photons. Diagonalizing Hamiltonian ^ in the 
basis \0n + I), \2n}} leads to the dressed states 

energies Q 



^o = 0, w± = i(A±n'), 

where 

n' = VA 2 +4g 2 (n+ 1) + 4ft 2 , 
and the dressed states 



(3) 
(4) 



|*o) = -sin#|ln) + cosi?|0n + l) 
|*_) = - sin y> [cos 0|1 n) + sin i?|0 ra+l)] 

+ cosip|2n) (5) 
|* + ) = cos y> [cos i? I Ira) +sin?9|0n + 1)] 

+ sint/?|2 n) . 



We have introduced the notation 
9 



sint? 



smip 



cos$ 



COS If 



Q' + A 



2fl> 



(6) 



(7) 



where g = g\Jn + 1/Cl. 

Initially, the cavity field mode is in a vacuum state 
and, therefore, we are especially interested in the time 
evolution of the 1 10) state, so we assume n = 0. Knowing 
the evolution for the dressed states ([5|), we can write 
the exact expression for this evolution of the initial state 
|10>(*) 



-iHt 



1 10) = |10)(t) = o(*)|10) + 6(t)j01) + c(t)|20), (8) 
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where 



a(t) = cos 2 t? 



+ /+(*)-§/-(*) 



b(t) 



sin i9 cos # 



A 



-! + /+(*)- 



(9) 



c(t) = *£/_(*), 



and 



(10) 



The perfect 7r pulse operation, defined by 1 10) — > |01), 
requires the condition g = 1, which means g = f2 and, 
therefore, we restrict the following considerations to this 
case only. Then, we have cos-d = sin d = l/y/2, and 
expressions (O take simpler form. 

It is evident from that the excited atomic state 
population shows oscillatory behavior. Since the qubit 
is coded into the lower states, the fidelity of quantum 
operations performed on the qubit is reduced when some 
amount of population is present in the intermediate level. 
To overcome this problem, the upper level is usually adi- 
abatically eliminated by choosing the detuning so large 
as to make the intermediate state population negligible. 
There is, however, an important drawback of this ap- 
proach — it is necessary to use very large values of the 
detuning to achieve fidelities of the quantum operations 
sufficient for quantum computation. It would require 
A/g > 640 to get the fidelity different from unity by 
amount smaller than 10 -5 . Such detunings lead to long 
operation times, which are very challenging for optical 
cavities. 

Here, we propose a scheme for fine tuning the quan- 
tum operations by using the fact that, in an ideal case, 
the population of the excited state evolves periodically in 
time, periodically approaching zero. By quantum opera- 
tion, we understand a unitary operation transforming a 
given initial state into another (desired) quantum state. 
The quantum operation is perfect if the desired state is 
produced with the fidelity equal to unity. The idea of 
fine tuning is to choose the operation time in a way as to 
complete the operation when the population of the inter- 
mediate state is zero. To achieve this goal, we choose the 
evolution time in such a way that f-(t = tk) = 0, which 
means 



(lu+ - uj-) t k = Q' t k = 2tt k, k = 1, 2, 3 ... , (11) 



and 



with 



/+(** 



exp 



i-(Q'-\A\)t k 



1 if A > 
-1 ifA<0 



(12) 



(13) 



Let us now require that, beside the relation (fTTj) . the 
evolution time satisfies the relation 

i(fi'-|A|)t M = |z, Z = l,2,3,...,2fc-1. (14) 

Both requirements for the evolution time can be satisfied 
provided the numbers k and I obey the relation 



Of 



2k 



I fi'-IA 



> 1, 2k > I, 



(15) 



which leads to the discrete set of detunings Ak,i given by 
the formula 



k.l 



2g 



2(f -1)+1 



(16) 



For large detunings 2k /I — 1 is large, and we can drop 
unity in the denominator getting simpler, but approxi- 
mate, relation. The discrete values of the detuning Ak,i, 
in units of 2g, are shown in Fig. [2] 




FIG. 2: Discrete values of jA^jj (in units of 2g) for initial 
values of k and I. 

Enforcing both condition (fTTj) and (fT4"]) gives us a dis- 
crete set of the evolution times for which solution ([5]) 
takes the form 

|10)(4,0 = \ [1 + (ie) 1 ] |10) (ie) 1 ] |01) . (17) 

The solution (|17[) is quite simple, and under appropriate 
choice of the numbers {k,l} it gives either the superpo- 
sition of the two initial states or one of the initial states. 
The important 7r operation we get for times 



for which 



t-n = tk,h f° r I even, 1/2 odd (18) 
e -iff^|l ) = _|01). (19) 



For the second basic operation, the 7r/2 pulse opera- 
tion, we should choose I odd, for which we have 



*7r/2 — tk,l, 



for I odd, 



(20) 
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and the solutions are 



IV. APPROXIMATE SOLUTIONS 



•|10) = 



i=*E|10> - ^|01), i=i odd 
i±ss|10> - ^|01), i=l even 



(21) 

Depending on the value of (Z — l)/2, we get from (|2ip 
one of the two orthogonal superposition states, which 
are maximally entangled (Bell states) of the atom-cavity 
system. 

The solutions presented above are a direct consequence 
of the periodic evolution of the system. The periodicity 
appears when the two frequencies u±, given by ©, are 
commensurate, i.e., their ratio is a ratio of integers. Tak- 
ing absolute values of the two frequencies we can distin- 
guish between the "fast" uj > and "slow" cu < frequency, 
where 



L0> = i(fi'±|A fc ,,| 



\w±\ ifA M >0 . . 



< 



From pd|) and (|14p it is easy to check that the frequencies 
lu < , uj > , and Q' = lo < + oj> are all commensurate, so the 
evolution is periodic. 

The period of the slow oscillation is equal to 



T = 



2tt 



4tt 



Q' - |A 



k.l\ 



tk,l 
I 



(23) 



and is related to the time tk,i for perfect quantum oper- 
ation. The period of oscillation of the intermediate state 
population is given by 



T' 



2n 

n 7 



2tt lo < 



T- 



l 

4k 



tk,l 
~k~ 



(24) 



The subsequent minima in the intermediate state pop- 
ulation are separated by T". If 4k /I is an integer the 
period of the system evolution is equal to T, and, simul- 
taneously, it is equal to an integer multiple (4k /I) of the 
period T' . When 4k /I is an irreducible fraction, the pe- 
riod of the system evolution is equal to 4kT' — IT, and if 
the fraction is reducible the period is reduced appropri- 
ately. 

For example, choosing k = 1 and I = 2, we have = 
and the resulting state is, according to (fil)|) , — 101), 
which is, up to the phase, an illustration of the perfect 
7r pulse operation. The operation is completed at time 
tk.i = kT' for which the population of the intermediate 
state is zero, so the fidelity of the operation is equal to 
one. Similarly, by choosing k = 1 and I = 1 we get, 
according to (f2"Tj) . the superposition state (1 +i)/2|01) — 
(1 — i)/2|10) which is generated with the fidelity equal 
to one. This is a perfect n/2 pulse operation, and the 
detuning in this case is equal to Ai ! i/(2^) = 0.8165. 

The two examples illustrate the idea of fine tuning of 
the quantum operations: choose one of the discrete values 
of the detuning Aj^j and corresponding operation time 
tk,i to complete the operation at a time when the popu- 
lation of the intermediate state is zero. 



So far we have discussed the ideal case, when there is 
no cavity decay, and the spontaneous emission from the 
atomic excited level is ignored. To make the system use- 
ful for quantum information processing it is necessary to 
have access to the quantum information stored in the sys- 
tem and the possibility of sending it over long distances. 
In the case of the quantum system under discussion it is 
possible, when one mirror of the cavity is partially trans- 
parent. Of course, transparency of the mirror leads to 
a damping of the cavity field mode. Moreover, sponta- 
neous emission introduces damping to the atomic system 
which spoils the desirable results of quantum operations. 
Unfortunately, in the presence of damping in the system, 
it is not possible to get exact analytical solutions; there- 
fore, we have to resort to some approximations. When 
the two decay rates k and 7 are small we can apply the 
perturbative methods to find the first order corrections 
to the dressed states energies as well as the state ampli- 
tudes. 

The first order corrections to the dressed states ener- 
gies lead to the following modifications: 

uj =u Q -iK , Q± = u± - i{n± + j±) , (25) 

where the damping rates kq, k±, and j± are given by 

«o = ~, K± = ^(lT<0, 7± = |(l±tf), (26) 

where d — A/57'. The solution ([5]) for the conditional 
evolution of the state |10), according to the first order 
perturbation theory with respect to both k and 7, is 
given, under the assumption Q = g (g = 1), by the fol- 
lowing formulas: 



a(t) 

b(t) 
c{t) 



X - {(1 + i 2 V d) e- K «* + (1 - i 2 V d) f + (t) 
- [d-t^t V (l + d 2 )] /_(t)} 
-{-e-^ + UV-id-i® f-(t)} 
^ {-* |e—* + < 2(1 + / +(t ) 



(27) 



where we have introduced the notation 9 = fl'/(2g), 
7 = 7/fi', k — n/ft', r) = R& 2 , and £ = 7/O 2 . Functions 
f± have the form (TIT))) except that the frequencies u>± are 
replaced by uj± from ([25]) . 

The solution (|27|) is valid as long as k and 7 are small. 
In fact, the real smallness parameters are 77 and £, so we 
require that both rj <C 1 and ^€1. In deriving l[2"7]) we 
also discarded terms proportional to the product ivy. On 
the other hand, the solution (|27|) is valid for any value of 
the detuning A. This means that it allows also for the 
resonant case of the ir pulse operation, when k = 1, 1 = 2 
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and A 1.2 = 0. The resonant case is usually not rec- 
ommended because of the spontaneous emission from 
the intermediate level. What is usually done to mini- 
mize the effect of spontaneous emission is the adiabatic 
elimination of the excited level. This requires large val- 
ues of the detuning A. It is seen from (|27|) that when 
1/0 = 2<7 /r2' <C 1, the amplitude c(t) of the excited level 
is small. Neglecting this amplitude by setting c(t) = is 
exactly what the adiabatic elimination is about. 

However, when the rate of spontaneous emission is 
small with respect to 2g, the approximate solution l|27p 
should describe properly the role of spontaneous emis- 
sion. This will be discussed later. 



A. Adiabatic elimination 

The standard procedure used to eliminate the influence 
of the exited atomic level on the result of quantum oper- 
ation is the adiabatic elimination of the excited level. It 
is realized by taking large detunings of the fields from the 
atomic transition frequencies, which makes the popula- 
tion of the excited state very small, and, in consequence, 
the state is assumed not to take part in the evolution. 
Assuming that 1/0 = 2g/Q! <C 1, which means that 
A/{2g) 3> 1 and the requirements for adiabatic elimina- 
tion are met, we can eliminate the state 1 20) from the 
evolution by putting in equations (|2"7|) c(t) = 0. We 
should also put 7 = (7 = 0, £ = 0) to get rid of the 7 
dependence in the other amplitudes. For large detuning 
we can also put |d| = 1. With all these substitutions we 
get the following solutions: 



o(t) = 



b(t) 



1 



- o [1 + 2ier) + (1 - 2ierf)e leu)<t ] 



— K(jt 



e 
1 

2 6 



-K0t a ^€UJ < t 



cos [ — a;<f I + 2i] sin I 2 W <^ 



/' sin [ 2^<^ 



(28) 



where now 77 = k|A|/(4<? 2 ). From (|28|) . one can calculate 
the time for the 7r pulse operation, which takes the form 



t v (K) = tk,l 



1 



2 1 

1 arctan — 

7r 2r\ 



and similarly for the 7r/2 pulse operation 

£ir/2(«) = tk,l 



, 4 1 

1 arctan 

7T 1 - 2r/ 



(29) 



(30) 



For large detuning co < = 2g 2 /\A\, 2r\ — k/uj < , and for- 
mulas (|2T)|) and ([3H)) are consistent with corresponding- 
formulas obtained when the adiabatic elimination is per- 
formed at the Hamiltonian level [sj. To be precise, the 
consistency is up to terms linear in k, because the solu- 
tions (|2"T|) . and consequently (j2"5)) , were obtained in the 



linear approximation with respect to k (or 77), i.e., under 
assumption r/ 1. As far as the linear approximation is 
valid, one can derive simplified formulas for the relative 
changes of the operation times, which take the form 

U(k) - t k ,i 877 

tk,i Ik 

^/2(«) ~ tk,i _ 4?7 

tk.l Ik 



(31) 



In the presence of damping the evolution is not unitary 
and the resulting state, under the condition that neither 
a photon from the cavity nor a spontaneous emission pho- 
ton are registered, is given by 



!*(*)} = 



!*(*)} 



(*(*)!*(*)} 



where 



!*!(*)> = 



-iHt 



1 10) 



(32) 



(33) 



is the unnormalized, conditional quantum state gener- 
ated with the effective Hamiltonian (JTJ) . The fidelity of 
the resulting state is equal to \(^(t)\lQ)(tk,i)\ 2 , where 
|10)(ifc,z) is a desired state given by (TTT)) . The so- 
lutions (|28|) are the amplitudes of the unnormalized 
state (1331) and must be normalized when we calculate 
state (1321) and the fidelity in the presence of damping. 
For simplicity we use the same notation for the normal- 
ized amplitudes a, b, c as for the unnormalized. 

It is well known that when damping is present in the 
system, the periodic behavior of the system is lost, and 
the maxima and minima are shifted. This is exactly what 
we observe here. The quantum states generated in the 
presence of damping are no longer the ideal states (fTO|) 
or (f2~Tj) . but if n/{2g) <C 1, the fidelity of the states gen- 
erated can be quite high. The best choice for th,i is 
that with the smallest possible values of I. An exam- 
ple is shown in Fig. [31 where the 7r pulse operation is 
illustrated. After adiabatic elimination of the interme- 
diate level the populations of the remaining levels oscil- 
late smoothly (broad lines), and the maximum is shifted 
from the time tk,i to However, without adiabatic 

elimination the populations are modulated with fast os- 
cillations, and there are fast oscillations of the population 
of the intermediate level seen at the bottom. With the 
parameters of Fig. [3jthe detuning A 3 i y2 /(2g) = 5.4321. 
If g/2ir = 16 MHz then A 31 , 2 /27r = 173.83 MHz, which 
is quite big, but the presence of the intermediate level is 
still visible, and it has significant influence on the fidelity 
of the resulting state. In the figures we illustrate the state 
evolution by plotting the state populations only, but the 
solutions are more general and give the state amplitudes 
that include also the phase information. The resulting 
state is a superposition of the basis states and the popu- 
lations do not fully characterize the superposition. More 
precise characteristic of the state created during the evo- 
lution is its fidelity. To make it clear we give examples of 
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CO 

a 



0.75 




0.25 



TABLE II: Normalized amplitudes of the quantum states for 
the 7r/2 pulse operation at various operation times (k — 31, 
I = 1, K /(2g) = 0.01) 



Time 


a 


b 


c 


Fidelity 


*fc,i 


0.5342+0.5351 i 


-0.4623+0.4632 i 


-0.0045+0.0029 i 


0.9948 




0.4584+0.5339 i 


-0.5371+0.4497 i 


0.0589-0. 1034 i 


0.9858 


tf 


0.4632+0.5407 i 


-0.5323+0.4579 i 


-0.0058+0.0029 i 


0.9999 



mation with respect to r\ it leads to the detuning 



FIG. 3: Comparison of the populations |a(t)| 2 , |6(t)| 2 , and 
\c(t)\ 2 obtained from Eqs. ([27]) and ([28]) (normalized) for k = 
31, I — 2, and n/(2g) — 0.01. Broad lines are populations 
after adiabatic elimination of the intermediate level. Vertical 
lines mark the times tk,i and t n (K). 



the state amplitudes and the fidelities in Tables HI and HT1 
For the example shown in Fig. [31 we present in Table U 
the amplitudes, calculated according to the approximate 
solution (|27|). for the state amplitudes at different oper- 
ation times. 



TABLE I: Normalized amplitudes of the quantum states for 
the 7r pulse operation at various operation times (k = 31, 
I = 2, k/(2q) = 0.01) 



Time 


a 


b 


c 


Fidelity 




-0.0029+0.1078 i 


-0.9940+0.0000 i 


-0.0162-0.0098 i 


0.9880 




0.0022+0.0066 i 


-0.9887-0.1012 i 


0.0372+0. 1043, i 


0.9877 




-0.0037+0.0064 i 


-0.9946-0.1015 i 


-0.0162-0.0133 i 


0.9995 



The time i/, which we denote as the fine tuning time, 
is defined as the time at which the closest to t ff (/t) max- 
imum of the population of the state 1 01} (or minimum 
of the population of the state 1 20) ) occurs. From insets 
in Fig. [3] it is clear that tf = tk,i + IT' , and this choice 
improves significantly the fidelity of the operation. Sim- 
ilar results are obtained for the 7r/2 pulse operation. An 
example is given in Table II, where the state amplitudes 
are presented for k — 31, / = 1, and n/{2g) — 0.01. In 
this case A 3 i i i/(2g) = 7.7784, and it is larger than in 
the previous case. This example is even more interesting 
because the fidelity at the operation time t v /2 ( K ) coming 
from the adiabatic elimination is smaller than the fidelity 
at time tk i for ideal case of periodic evolution. The fine 
tuning, which this time gives tf — tk,i + 3T', again im- 
proves the fidelity. From Fig. [3] it is clear that the time 
£„•(/«) is different from the time tk,i for the ideal case. 
This difference suggests another possibility of fine tuning 
of the quantum operations, which consists in taking the 
operation time calculated from (|29[) as a new tk.i and ad- 
just the detuning appropriately to get a maximum of the 
fast oscillations at the new tk,i- Within linear approxi- 



|Afc,j(«OI 



|A fc 



k.i 



4?7 . 2 
— sin 

47T 



(34) 



By adjusting the detuning according to (j34j) we get the 
situation illustrated in Fig. 0J Now, the fidelity at time 




0.5 1 

FIG. 4: Same as Fig. [3] but with A k ,i( K ) S iven b y <G2J- Ver- 
tical line marks the ^(/t). 



t n (n) takes the value 0.9992. The fine tuning of the quan- 
tum operations can be performed in both ways: by find- 
ing the best operation time and/or adjusting the detun- 
ing. The values given by (fTtj]) for the ideal case are 
a good starting point. 

Examples given here show clearly that even for quite 
large values of detuning, the presence of the intermedi- 
ate state is still visible, but the fidelity of the result- 
ing states can be significantly enhanced when the opera- 
tion time and/or the detuning are adjusted appropriately. 
This, however, requires knowledge of the solutions for the 
three-level system, and cannot be done when the adia- 
batic elimination of the intermediate state has already 
been accomplished. 

Making the adiabatic elimination on the Hamiltonian 
level, consisting in diagonalizing the 2x2 resulting Hamil- 
tonian without resorting to perturbation theory, leads to 
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the following results 



a(t) 



b{t) 



c;2 



(eu <+ire)t 



COS ( — CL><(/c)t 



H 7-7 sin ( ^w<(«)t 

" ; — — — ^'111 I -w<(«;)f 
w<(k) V 2 



,(35) 



where cj<(k) = y c<j 2 — k 2 is the modified by k frequency 
w< . Since the modification is at least quadratic in k, by 
keeping only the linear terms in (|35j) , we easily reproduce 
formulas (|2"8]l. When corrections due to k 2 become im- 
portant, the populations of the two states oscillate with 
the modified frequency lo < (k) instead of w<. The opera- 
tion times are then given by 



tn{K) = 

r 



2 1 2 
1 + - 1 arctan 

I \ 7T 



1 Tt/2 



tk,l 



r 



1 / 4 
1 1 arctan — 

l V t 1- VT^T 2 



where r is the ratio given by 



r = uj < (k)/uj < = y/l - (k/w<) 2 



(36) 



(37) 



Since u> < — 2g 2 /\A\ decreases as |A| increases, the ra- 
tio r for a given value of n can significantly differ from 
unity. This is true for very large detunings, i.e., for large 
values of k. It is important, however, to have k,/lj < < 1, 
which puts some restrictions on the possible values of the 
detuning. 



B. Small detunings and spontaneous emission 

In the case of small detunings, adiabatic elimination of 
the intermediate state is not possible, but solutions (|27|) 
are valid for arbitrary detuning, that is, also for small 
detunings and resonant cases. In order to check how 
close to unity the fidelity of the n pulse operation can 
be, we have investigated this problem for the parame- 
ters (n,g,n)/2n = (16,16,0.05) MHz. This means that 
n/{2g) = 0.0015. Of course, for nonzero cavity damping 
k, we should choose the smallest possible values for I as to 
complete the operation in the shortest possible time. For 
example, for k = 1, I = 2 we get from (|T6)) Ajy = and 
we obtain, after time tjy, the state 1 01) with the fidelity 
different from unity by less than 10~ 6 . This is the reso- 
nant case, which is usually not recommended because of 
the spontaneous emission from the excited atomic state. 

Similarly for the tt/2 pulse operation, for k = 1, I = 1 
we get A/(2g) = 0.8165, which means that for g/2TT = 
16 MHz, A/2rr = 26.128 MHz, and for such a value of 
A we can reproduce the state (|2"Tj) with the fidelity that 
differs from unity by less than 10 -6 . 

Choosing k = 4 and I = 2, we get the ir operation with 
A ky i/(2g) = 1.6036 (51.314 MHz for 3 /2tt = 16 MHz), 



and we have a noticeable value of detuning. For k/2tt = 
0.05 MHz, we still get fidelity 0.99997. 

If spontaneous emission is ignored, one can expect very 
good performance of quantum operations at small detun- 
ings. The approximate formulas (|2T[) take into account 
spontaneous emission and they are valid whenever the 
spontaneous emission rate 7 is small with respect to fi'. 
We can thus use equations (|2"T)) to calculate the state am- 
plitudes in the presence of spontaneous emission. Assum- 
ing the parameter values (f2, g, k, 7)/2tt = (16, 16, 0.05, 1) 
MHz, we have n/(2g) = 0.0015 and 7/(2.9) = 0.03 and, 
for k = 1 and 1 = 2, get the situation illustrated in Fig.[5j 
where the state populations are shown. With these val- 



3 



0.75 



0.5 



0.25 



\ / 


\ " /~ 

\ / 


\ / 


\ / 


\ / 


\ I 


\ / 




\ / 
\ / 


\ / 


\ / 
\ / 


\ / 
\ / 
\ / 


\ -- ' 
\ v 


\ / 




















: / \ 




/ / X 



0.5 1 

t/tk,l 



1.5 



FIG. 5: Evolution of the atom-field states populations (nor- 
malized): |a(t)| 2 (solid line), \b(t)\ 2 (dashed line), |c(i)| 2 (dot- 
ted line) for k = 1, I = 2, ft/(2p) = 0.0015, and -y/(2g) = 0.03. 
Time is scaled in units of tk,i- 

ues of parameters the state 1 1 ) is generated with the 
fidelity still equal to 0.9989, which is quite good. For the 
7r/2 pulse operation the situation is even better and the 
fidelity takes the value 0.9993. Of course, spontaneous 
emission rapidly deteriorates the quality of states gener- 
ated at small detunings, but it is still possible to obtain 
remarkable fidelities. 



V. NUMERICAL RESULTS 

Approximate analytical results presented in previous 
section are valid only for sufficiently small values of k 
and 7. They illustrate the idea of fine tuning of quan- 
tum operations, but looking, for example, at Fig. [4] it is 
clear that it should be feasible to get even better results 
when more precise tuning is carried out. To this end, nu- 
merical calculations are necessary. The effective Hamil- 
tonian |T]) can be diagonalized numerically and we can 
find nonunitary evolution governed by this Hamiltonian. 
The numerical solution allows for adjustment of the de- 
tuning and corresponding time allowing for further im- 
provement in the fidelity For example, for the situation 
illustrated in Fig. [H numerical optimization gives the fi- 
delity 0.9997 if we take A 3 i j2 (k)/(2 5 ) = 5.2409 instead of 
A 31 , 2 (K)/(2g) = 5.2380 given by ® and A 31 , 2 /(2 ff ) = 
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5.4321 given by (fTB"]) . which for g/2n — 16 MHz gives 
for the dctvmings A 31i2 /27r values 167.71, 167.62, and 
173.83 MHz, respectively. This example shows that for 
n/(2g) — 0.01 and 7 = 0, analytical formulas give quite 
accurate results, and precise numerical tuning does im- 
prove the results but only slightly. 

The situation changes dramatically when k and 7 are 
not so small as to justify linear approximation made in 
derivation of the analytical formulas. In such a situation 
the only reliable solutions are numerical solutions. To vi- 
sualize the role of the two damping rates in the evolution 
we plot in Fig. [6] the populations of the quantum states 
for the value of k/ (2g) = 0.05 and 7 = 0. The time t v (K) 
is now calculated according to ([55]) . The amplitudes of 
the fast oscillation are quite big, and the ir pulse oper- 
ation is far from being perfect. The detuning A 31 2 is 
numerically tuned to fit the maximum of population and 
takes the value A 31t2 /(2g) = 4.4491. We find that the 



value 7/(2*7) = 0.03 for the spontaneous emission rate, 
which gives for g/2ir — 16 MHz the value 7/27r = 1 MHz. 
It turns out that the presence of spontaneous emission 
has a stabilizing effect on the evolution. The fast oscil- 
lations are damped significantly, and, surprisingly, the 
fidelity of the state generated becomes higher. For the 
parameters of Fig. [7] we find the fidelity 0.9994. The rea- 
son for such behavior can be understood from Eq. (|26p . 
which shows the difference between k± and 7± when the 
detuning is large and \d\ — *■ 1. The two damping rates act 
in a sense in opposite directions: when fast oscillations 
are weakly damped by k they are strongly damped by 7, 
and vice versa. The effect is more pronounced for higher 
values of the detuning, i.e., for higher values of k. Sim- 
ilar effects can be observed for the ir/2 pulse operation. 
These results show that in some situations spontaneous 
emission can play a positive role in performance of quan- 
tum operations. 
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FIG. 6: Numerical solutions for the state populations 
for K,/(2g) — 0.05, "//(2g) — 0, and numerically tuned 
A 31 , 2 /(2<?) = 4.4491. 

fidelity in this case is equal to 0.9660 and is, of course, 
much worse than it was for k/ (2g) — 0.01 shown in Fig. [4] 
In both cases the spontaneous emission is not taken into 
account. In Fig. [7] we present the same situation but 
when the spontaneous emission is included. We take the 
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FIG. 7: Same as Fig. ©but for 7/(25) = 0.03 



A. Nonrectangular pulses 

In real experiments the field cannot be switched on 
and off abruptly, so the rectangular pulses are rather not 
realistic. There are always a finite rise time and a finite 
fall time of the pulse. The pulse has definite shape, and 
usually what we observe in experiments depends on the 
pulse shape. Here, we address the effects of the pulse 
shape on the fine tuning discussed above. 

With the Hamiltonian fl}, assuming that — > gF(t) 
and g — > gF(t), where the function F(t) describes the 
pulse shape, we get from the Schrodinger equation (units 
scaled to 2g) the following set of equations for the state 
amplitudes 



d(t) 



-F(t)c(t), 



(38) 



On resonance, without damping, the system [[38]) has the 
solutions 



b(t) = -nb(t)--F(t)c(t), 
m = ~\F{t)a{t)- l -F{t)b{t)-{ 1 + l A)c{t). 



a(t) 
b(t) 
c(t) 



"[! + /+(*)] , 
"[-! + /+(*)] 

T2 f - {t) - 



(39) 



where 



/±(*) 



2^ XP 



± exp 



-7lL Hf)it ' 



(40) 
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The solutions (|39[) depend only on the pulse area 
L F(t')dt' , so for any shape F(t), with the same area, 
the resulting state after time t is the same. 

We do not know analytical solutions for the nonreso- 
nant case, even without damping, but the system ([38]) 
can be easily solved numerically for a given pulse shape 
F(t). The solutions lack their periodic character and 
the quantum operation time depends on the pulse shape. 
The situation is thus more complex than for the rectan- 
gular pulses discussed in previous sections. Nevertheless, 
even in this case the idea of fine tuning is useful for find- 
ing the operation time. We illustrate the problem with 
few examples. Let us consider two different pulse shapes: 
the trapezium shape and the sine square shape, which are 
defined by the functions: 



F{t) 



for < t < t r 
for t r < t < t 



p-tf 



s^- for tp-tf <t <t p 



where 



1 



tr+tf 



(41) 



(42) 



is the scaling factor which makes the pulse area for the 
pulse duration t p to be the same as the area of the rect- 
angular pulse of the same duration and unit height. The 
times t r and tt are the rise and fall time of the pulse, 
respectively. The normalization to the same pulse area 
ensures, on resonance, the same final state at time t p 
after the pulse for different pulse shapes. 

The sine square pulse shape is given by the function 
(normalized to the same area) 



F(t) 




for < t < t p 
otherwise 



(43) 



The two pulse shapes are of different character. The 
trapezium is close to the rectangle when the rise and fall 
times are short. The sine square pulse is much narrower 
and changes smoothly over the whole duration time t p . 
For the rectangle pulse we have F(t) = 1. 

To find approximately the proper operation times for 
pulse excitation, we use the following trick. We replace 
the oscillation frequency f2' of the intermediate state pop- 
ulation with an "average frequency" 



n' = VA2 + 2 -»• Q' = J A 2 + 2F 2 (t) , (44) 



where 



1 



which gives 



F^t) 



F*{t) = ^l F 2 (t)dt 
1 p Jo 



s(4 — s)/3 for trapezium 
3/2 for sine square 



(45) 



(46) 



with s given by (|42J>. Applying such a replacement 
we can calculate the "period" T" using l|24p and the 
time tk,i for the considered operation using Aj^j defined 
by (|16p . The pulse is characterized by one parameter — 
the mean square amplitude defined by (|45|) . For pulses 
with smooth rise and fall times, the evolution is no longer 
periodic, so the time X" is not really the period, but for 
small values of I (I — 1,2), the time tk,i calculated in this 
way is still a good approximation for the optimal opera- 
tion time. It can at least be treated as a good starting 
point for further optimization. 

To illustrate this situation we plot the evolution of the 
final states populations for two different pulse shapes. In 
Fig. Owe present the evolution of final state populations 
for the trapezium shape with t r /t p = tf jt p = 0.1 and the 
sine square pulse of the same duration. The other param- 
eters are chosen as k = 3, I — 2, n = 0, 7 = 0, and the 
pulse duration t p = tk,i with tkj calculated with Q' ad- 
justed according to (144j) . The small value of k — 3 gives 
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FIG. 8: Populations of the resulting states for pulse excitation 
with the trapezium pulse (upper figure) and sine square pulse 
(lower figure) for k = 3, I = 2, k — 0, and 7 = 0. The meaning 
of the lines is the same as in Fig. [5] 

the detuning A^^/i^g) = 1-2649 which is rather small, 
and in this case, as is seen from Fig. [8j the desired state 
is not of good quality. The fidelity, assuming the pulse 
duration t p = tk,i, is equal to 0.8176 for the trapezium 
pulse and to 0.9529 for the sine square pulse. This means 
that the approximation for t^j is not very good for small 
detunings. Small detunings mean the situation is close to 
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resonance, and the solutions should be close to periodic 
solutions ([551) . Actually, the oscillations are still visible 
in the figure. The fidelities of the generated states can be 
improved by adjusting the pulse duration. If the pulse 
duration is increased by a factor of 1.103, for the trapez- 
ium pulse, we get the fidelity equal to 0.9681, and this is 
the best result for the trapezium with A3. 2. In the case 
of the sine square pulse, by increasing the pulse length 
by a factor of 1.087, we get a fidelity better than 0.99999. 
In this respect, for small detunings, the sine square pulse 
appears to be better than the trapezium pulse. We have 
to remember, however, that the results refer to the ideal 
situation, where there is no damping. 

In Fig. [9] we illustrate the situation for k = 31 and 
1 = 2. In this case we have detuning A^i^/i^g) = 5.4321, 
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FIG. 9: Same as Fig. but for k = 31 and I = 2. 

and the pulse duration time t p = tk t i calculated according 
to the above procedure fits much better for producing 
desired state. The fidelities obtained with such t p are: 

0. 99994 for the trapezium shape and 0.99899 for the sine 
square pulse. The results are remarkable and confirm 
that the procedure works pretty well for larger detunings, 

1. e., higher values of k. The result for the trapezium pulse 
is optimal for the detuning Aai^, while the result for the 
sine square pulse can be still improved when the pulse 
duration t]~,i is multiplied by a factor of 1.02 giving a 
fidelity better than 0.99999. Since the evolution is not 
periodic as it is in the case of rectangular pulses, the 



procedure will not work so well for higher values of /. The 
results obtained give evidence that for pulse excitation 
the idea of fine tuning of the quantum operations is also 
applicable. 

In a real experimental situation, we have to take into 
account the damping rates ft and 7. In analogy to the 
results for rectangular pulses, we can expect that the 
pulse duration should be increased to get desired states 
with good quality. The approximation described here 
does not give any hint as to how long the pulse should 
be, but starting from the ideal case it is easy to adjust 
the time numerically. For example, for the sine square 
pulse, taking n/(2g) = 0.05, 7/(2g) = 0.03, and other 
parameters as in Fig. [51 we get a fidelity better than 
0.99999 for the time 1.154 Even in this case the 
approximate value of the pulse duration t p = tk,i is a 
good starting point for numerical optimization. 



VI. CONCLUSION 

In conclusion, we have studied the effect of the popu- 
lation of the intermediate state on two basic operations 
performed via the two-photon Raman transition in the 
atomic A system. Full dynamics of the three-level sys- 
tem is taken into account without adiabatic elimination 
of the intermediate level. The dynamics studied are con- 
ditioned on no photon detection emitted by spontaneous 
emission and no photon detection leaking out of the cav- 
ity. We have found a discrete set of detunings leading 
to perfect tt pulse and tt/2 pulse operations. It is shown 
how to make the population oscillations for improving 
the quantum operation precision useful. In the ideal case 
of « = and 7 = 0, the population of the intermediate 
state 1 2) oscillates, becoming periodically zero and giv- 
ing perfect operations. For nonzero, but small values of ft 
and 7, the fidelity can achieve values very close to unity. 
Approximate analytical solutions obtained in first order 
perturbation theory are given, including both the cavity 
damping rate k as well as the spontaneous emission rate 
7. Numerical results show that at certain circumstances, 
the spontaneous emission can play a positive, stabilizing 
role and improve the performance of quantum operations. 
The idea of fine tuning can also be applied for the pulse 
excitation with pulse shapes different than rectangle. 

For realistic, moderate values of the detuning, the in- 
termediate state has significant influence on the reliabil- 
ity of the quantum operation, but its destructive role can 
be limited by the fine tuning presented here. The pos- 
sibility of using small or moderate values of detunings, 
without reducing the operation fidelity by the nonzero 
population of the excited state, is also important because 
the operation time is a function of the detuning, and 
short operation times are preferred to minimize effects of 
dissipation. 
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